Gaussian Tunneling Model of c-Axis Twist Josephson Junctions 
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We calculate the critical current density j/ for c-axis Josephson tunneling between identical high 
temperature superconductors twisted an angle (j>o about the c-axis. We model the tunneling matrix 
element squared as a Gaussian in the change of wavevector q parallel to the junction, {|f(q)p) oc 
exp(— q^a^/27r^(T^). The j/(<^o)/J'/(0) obtained for the s- and extended-s-wave order parameters 
(OP's) are consistent with the Bi2Sr2CaCu2 08+i data of Li et al, but only for strongly incoherent 
tunneling, > 0.25. A dj.2_j^2-wave OP is always inconsistent with the data. In addition, we show 
that the apparent conventional sum rule violation observed by Basov et al. might be understandable 
in terms of incoherent c-axis tunneling, provided that the OP is not d^2_„2-wave. 
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Recently, there has been a resurgence of interest 
in the symmetry of the order parameter (OP) in the 
high temperature superconductors (HTSC). Al- 
though many phase-sensitive experiments were inter- 
preted as giving evidence for an OP in YBa2Cu307_a 
(YBCO) that was consistent with the ci3.2_j,2-wave 
form, the same type of phase-sensitive experiments 
on the electron-doped HTSC Ndi.85Ceo.i5Cu04_y and 



(PCCO) also were interpreted in 
in apparent contra- 
suggesting a possible 



Pri.85Ceo.i5Cu04 
terms of a d^2_y2-wave OP, 
diction with other results 
problem with the former experiments. In addition, 
new phase-sensitive experiments on Bi2Sr2CaCu208-i-5 
(Bi2212) gave strong evidence that the OP contains an 
isotropic component, since it has a non- vanishing aver- 
age over its Fermi surface. First, c-axis Josephson 
tunneling between Bi2212 and either Pb or Nb demon- 
strated that an isotropic s-wave component of the OP 
exists in Bi2212. j^J^ Second, c-axis tunneling across the 
junctions of Bi2212 intercalated with HgBr2 has been 
studied using mesas. |^ In these experiments, increas- 
ing the c-axis spacing by 6. 3 A increased the normal state 
resistivity by a factor of 200. In the superconducting 
state, Rn and Ic changed by comparable factors, but their 
product IcRn ~ 10 mV, about half the optimal value ex- 
pected in the Ambegaokar-Baratoff (AB) model of purely 
incoherent c-axis tunneling between identical s-wave su- 
perconductors. HJl^ Such behavior is very difficult to 
understand in terms of a d2.2_y2-wave OP. []lT|] 

More directly, a new phase-sensitive experiment which 
can test the symmetry of the OP over the entire range 
of temperatures T below the superconducting transition 
temperature Tc was performed. |^ In this experiment, a 
single crystal of Bi2212 was cleaved in the a6-plane, the 
two cleaves were twisted a chosen angle 4>o about the c- 
axis with respect to each other, and fused back together. 
After lead attachment, the critical currents //(T) and 
/^(T) across the twist junction and single crystals were 
measured, and their ratio was found to be directly pro- 



portional to the ratio A'^ /A^ of their areas at T/Tc = 0.9. 
1^ Those authors then claimed that (a) the intrinsic junc- 
tions and the twist junction behaved identically, (b) the 
c-axis tunneling is strongly incoherent, and (c) the OP 
contains an isotropic component, but not any purported 
(i2.2_j^2-wave component for T < Tc, except possibly be- 
low a second, unobserved phase transition. |^ 

Since then, the group theoretic arguments upon which 
conclusion (c) were based have been published. ||l2| In 
addition, an exact calculation of the possible roles of co- 
herent c-axis tunneling has been presented. [|3| For the 
tight-binding Fermi surface generally thought to be appli- 
cable to Bi2212, |ljjl^, it was shown that such coherent 
tunneling was inconsistent with the data, |^ even for an 
isotropic s-wave OP. Since the claim (a) of Li et al. is 
just a statement of their experimental observations, and 
the claim (b) is clearly correct in the limit of purely in- 
coherent tunneling, it remains to quantify precisely just 
how incoherent the tunneling must be in order to fit the 
data. 

Since Bi2212 behaves as a stack of weakly-coupled 
Josephson junctions, p6[ | the critical current density Jc 
across each junction may be evaluated by neglecting the 
couplings between the other junctions. Specifically, 
j/ across the twist junction between layers 1 and 2 is 
given by 



J/ = |4eT^(/'^(k-k')Fi(k^)Ft(k/))|, 



(1) 



where oj represents the Matsubara frequencies, /''(k— k') 
is the spatial average of the tunneling matrix element 
squared, (. . .) represents two-dimensional integrals over 
each of the two first Brillouin zones (BZ's), and the 
wavevectors kj and kj are obtained from k = {kx,ky) 
and k' — {k'^,k'y) by rotations of ±</)o/2 about the 
c-axis, respectively. The anomalous Green's functions 
F„ = A„/[o.2 + + |A„|2] and = F*, where A„(k) 
and ^„ (k) are the OP and quasiparticle dispersion on the 
nth layer, respectively. For Bi2212, we assume the quasi- 



1 



particle dispersions ^„ have identical tight-binding forms, 
differing only by the wavevector rotations, 



^„(k) = —t[cos{kxa) + cos{kya)] 

+t' cos(fcj,a) cos{kya) — /i. 



(2) 



where we take t =306 meV, t'/t = 0.90 and fi/t = -0.675 
to give a good fit to the Fermi surface of Bi2212, for which 
^„(ki?) = 0. These values are slightly different from those 
used previously, [jl^ A plot of this Fermi surface is shown 
in Fig. 1. We study three OP's. These are the isotropic s- 
wave, a constant, the dx2_y2-wave, topologically equiva- 
lent to cos(fca;a) — cos(A:j,a), or an extended- s-wave, which 
we take to be the absolute magnitude of the particular 
dx2-y2-wa,ye form. To obtain the particular d2.2_j,2-wave 



form we used the repulsive interaction of the form [|l8yi9| 



(3) 



Q=(±l,±l)-!T/a 



where Vq — 556 meV and F — 0.1. The BCS equation is 
solved for A(k) in terms of V^(k — k'). The extended-s- 
wave OP was taken to have the absolute magnitude of the 



da;2_„2-wave form obtained from Eq. 



The isotropic 



s-wave OP was taken to have the maximum magnitude 
of it. We included the umklapp terms which occur with 
rotated Fermi surfaces, but the effects of doing so were 
very small in all cases we studied. 




FIG. 1. Schematic plot of the Bi2212 Fermi surface used in 
these calculations. The concentric circles represent the por- 
tions of the BZ for which the tunneling strength is within 1/e 
of its maximum value, for dimensionless Gaussian parameters 
of magnitude 0.005 and 0.05, as indicated. 

Previously, we investigated the case in which f'^(k — 
k') (X (5(k — k'), appropriate for coherent tunneling. jT^ ] 
We found that it made very little difference whether or 
not one included umklapp terms in the integrals in Eq. 
(1). However, we found that for Fermi surfaces similar 
to that pictured in Fig. 1, it was not possible to obtain 
a quantitative fit to the data. We then argued that 

for extremely incoherent tunneling, with /''(k— k') = /g', 
a constant, the s-wave and extended-s-wave OP's could 



both fit the data equally well. However, we did not dis- 
cuss the intermediate case in any detail, and that is our 
purpose here. We therefore take the tunneling matrix ele- 
ment squared to be a Gaussian in the momentum change 
parallel to the junction. 



f'{k - k') = f,' exp[-(k - k'fay27T'a% 



(4) 



Since cr = and cr — s- cx) result in purely coherent and in- 
coherent tunneling, cr is a dimensionless parameter quan- 
tifying the incoherence of the tunneling. In Fig. 1, the 
shaded concentric circles for cr^ = 0.005 and 0.050 are 
pictured for the case in which the tunneling on one side 
of the junction occurs from the intersection of the Fermi 
surface and the hue connecting the (0, 0) and (tt, tt) points 
in the BZ. In the shaded regions, |k — k'| < any/2/a, so 

that 1/e </''(k-k')//6'<l. 

In the case of a circular Fermi surface cross-section, as 
might be expected for electron doped HTSC, we could 
write k = fc_F(cos (/)k;, sin (/)k), etc. and Eq. (^) can be 
rewritten as 



/■^('/>k,0k') = /"'(7)exp[7cos(0k - M], 



(5) 



where 7 = {akp/ira)'^ and /''{"f) = /o'exp(— 7). This 
form is identical to that of Graf et al, provided that 
/q — exp(7)/[Tx/o(7)]j where In{z) is a Bessel func- 
tion and 1/tj_ is an effective interlayer tunneling rate. 
The denominator contains the normalization factor /o(7)- 
Then just below T^, A{T) 0, and 



Jcc = C^ol(/''(0k,0k')Aic(0kJA2c(0k'J)0.,0„|, 



(6) 



where Cq — eni^/{4Tc), m is the in-plane effective mass, 
C — s,d,e indexes the OP's, and (/)k+ = 4>\l + 0o/2, 
4>-k'_ = 0k' - 0o/2. These OP's are A„s = Aq, A„d((/)) = 
Aocos(20), and A„e((/)) = |Aocos(20)|, respectively. We 
find 



•Jed 

and 
where 

and 



CoA2|cos(20o)|/d(7)/T±, 
CoA2/,(7,0o)/t_i, 

^2(7) 
2/0(7) 



rt=0 



- (5n,o)cOs(4n(^o)/4n(7) 

7r2(4n2- 1)2/0(7) 



(7) 
(8) 

(9) 
(10) 

(11) 



In this simple model, Jce/Jcs — /e(7, 0o) at Tc, 
which is periodic in (po with period 7r/2, satisfying 
/(7,±7r/2— 4>o) = f{'j,4'o)- In the incoherent limit, 
/e(0, (po) = (2/7r)2, a constant. However, in the coherent 
limit /e((X3, 0o) = [| sin(20o)| + (^/2 - 2|0o|) cos(20o)]/^ 
in the domain |(/)o| < tt/2. /e(cxD,0o) has maxima of 1/2 
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at (j)o = n7r/2 and minima of I/tt at (po = {2n + l)7r/4 
for integer n. The 4>q and 7 dependencies of fei"f, 0o) are 
shown in Figs. 2 and 3. 
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FIG. 2. Plots of /e(7,</)o) = rx Jci((/)o)/CoAg just below n 
for the extended- s- wave OP Ao| cos(2<j!>k)| on a circular Fermi 
surface, as a function of (po in degrees, at the 7 values 10"^, 
where a; = 0, 1/2, 1, 3/2, 2, and for the coherent limit 7 00. 

However, for all 7 values, J;;^ oc | cos(2(/)o)|, which van- 
ishes at 00 = 45°, as shown in Fig. 4. In addition, 
J^(0o = 0) differs from the s-wave case by the addi- 
tional factor fdi'y) = -^2 (7)/ [2/0 (7)]. In the coherent 
limit, 7 — > 00, fdin) — * 1/2, but in the incoherent limit 
7 <C 1, fd{l) — ^ 7^/16. /e(7, 0) also decreases from 1/2 
at 7 — > cx) with decreasing 7, but not nearly as dramat- 
ically, approaching the constant value (2/7r)^ as 7 0. 
The functions fd{l) and /e(7, 0) are compared in Fig. 3, 
where they are plotted as functions of 1 for clarity. 
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FIG. 3. Plots of T^Ji /Co^l at and 0o = 0, which is 
/e(7, 0) and fd{'y) for the extended-s- and d-wave OP's on the 
circular Fermi surface, respectively, versus I/7. 

We now consider the full calculations with the Fermi 
surface shown in Fig. 1, and Gaussian tunneling, Eq. (jj). 
We performed explicit calculations at T/Tc = 0.5, 0.9. 
In Figs. 4-6, we plotted our j/((/)o)/J/(0) results at 
T/Tc = 0.5 for the d^2_y2-'wave, extended- s-wave, and 



isotropic s-wave OP's, respectively, for 0° < (/)o < 45°. 
Results for T/Tc = 0.9 are very similar. All results are 
periodic in (f>o with period 7r/2 and symmetric about 
4>o = ±7r/4. Note that this presentation is equivalent 
to j/((/)o)/Jf, since is the critical current density 
across the untwisted single crystal, which is itself a stack 
of essentially equivalent Josephson junctions. Hence, 
^ci^) — Jc ■ III sach case, we show the results for 
cr2 = 0, 0.005, 0.05, 0.25, and 0.50. 
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FIG. 4. Plots of j/((^o)/J/(0) at T/Tc = 0.5 for twist 
junctions with the Fermi surface shown in Fig. 1, and a 
d^2_y2-wsive OP. The dimensionless Gaussian parameters 
are indicated. The curves for = 0.25 and 0.50 are nearly 
indistinguishable from the function | cos(2<^o)|. 
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FIG. 5. Plot of {(f>o) / Jc (0) at T/n = 0.5 for twist 
junctions with the Fermi surface shown in Fig. 1, and an 
extended-s-wave OP. The dimensionless Gaussian parameters 
are indicated. 

For each OP, coherent tunneling causes j/(0o)/>^/(O) 
to decrease sharply as (j>o increases from 0°. This is be- 
cause the tight-binding Fermi surfaces are only identi- 
cal for (j)Q = 0°, and the overlap of the Fermi surfaces 
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changes dramatically from a continuous curve to a set 
of four points for > 0°. In addition, there are local 
maxima in J/((/)o)/'^/(0) at about c^o ~ 20°, which are 
absent for a circular Fermi surface cross-section. For this 
twist angle, the rotated Fermi surfaces begin to overlap 
at eight points. This increased overlap is responsible for 
the more gradual decrease of J/(</'o) at larger val- 
ues. We also note that the coherent tunneling curves 
for the extended-s- and d-wave OP's are nearly equal for 
00 < 30°. This is also true for the — 0.005 curves. 
However, the d-wave curves are distinctly different from 
the s-wave and extended-s-wave curves in the vicinity of 
(po = 45°, where the d-wave curves all vanish by sym- 
metry. In addition, the d-wave curves for = 0.25 and 
a'^ = 0.50 are essentially identical and quantitatively in 
agreement with the function | cos(2(/<o)|- This is the same 
00 dependence that we obtained analytically for all 7 in 
the case of a circular Fermi surface cross-section. 
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IcRn of the critical current times the resistance across a 
junction from the AB limit {IcRn{T)\AB) for the case of 
an extended-s-wave OP as cr — s- cx). We first define 




FIG. 6. Plot of Jc iM / Jc (0) at T/Tc = 0.5 for twist junc- 
tions with the Fermi surface shown in Fig. 1, and an isotropic 
s-wave order parameter. The dimensionless Gaussian param- 
eters are indicated. 

In comparing the cases of the isotropic and extended- 
s-wave OP's, we see that the main differences appear 
for rather coherent tunneling, for which the extended-s- 
wave curves have a more dramatic 0o-dependence. How- 
ever, both OP's are consistent with the data of Li et al. 
for cr^ > 0.25, and neither one is consistent with the 
data for cr^ < 0.05. By examining Fig. 1, we see that 
= 0.25 corresponds to a tunneling strength satisfying 
1/e < /"'(k-k')//6' < 1 over 39% of the first BZ, which 
is strongly incoherent. Thus, these calculations imply 
that it is impossible to distinguish the isotropic-s-wave 
from the extended-s-wave OP's in the c-axis twist junc- 
tion experiments. We can only say that there must be an 
isotropic s-wave OP component for T < Tc, which is pro- 
jected out from the OP anisotropy, if any, by averaging 
over a large portion of the BZ. 

Finally, we have studied the reductions in the product 



IcRniT) - C(r/T,)/ei?„(T)| 



AB, 



(12) 



where the AB limit curve corresponds to the isotropic 
s-wave case with cr —s- 00 for any Fermi surface. A few 
limits can be investigated analytically. For a circular 
Fermi surface cross-section, we can analytically evalu- 
ate IcRn for this model both at Tc and as T ^ 0. We 
find C(l) = (2/7r)2 w 0.405, as above, and C(0) = 
(2/7r)2ln4« 0.562. 

Similar results can be obtained for different Fermi 
surfaces. For the Fermi surface shown in Fig. 1, we 
find numerically that C(0.9) = 0.416, C(0.5) = 0.465, 
and C(0) = 0.572. For the shghtly different Fermi 
surface studied elsewhere, ||r^ , |l7|| with t'/t — 1.3 and 
fi/t = -0.6, we found C(0.9) = 0.400, C(0.5) = 0.450, 
and C(0) = 0.578. Thus, the result of Yurgens et al. 
that IcRn{T = 0) w 10 mV for the HgBrs-intercalated 
Bi2212 is in rather good agreement with that expected 
for an extended-s OP. Q 

We remark that there have been a number theories of 
the HTSC that relied upon coherent tunneling in the c- 
axis, with a k-dependent matrix element, between each 
of the layers. |^,^ This model was shown previously 
to give the worst agreement with the experiment of Li et 
al. of all of the coherent tunneling models studied, even 
for an s-wave OP. [|6|p^ Here we showed that the exper- 
iment is easy to understand if the interlayer tunneling 
is strongly incoherent, provided that the OP integrates 
to a non-vanishing value on the two-dimensional Fermi 
"surface" . 

Although there have not yet been many infrared re- 
fiectance measurements on Bi2212, the available exper- 
iments strongly suggest that the c-axis tunneling is not 
metallic. j2^ In that work, not only was no Drude edge 
above 30cm^^ seen in the c-axis conduction in the nor- 
mal state, none was seen in the superconducting state as 
well. In addition, there is now strong evidence at all un- 
derdoped HTSC have incoherent c-axis normal state tun- 
neling, with the sole exception of YBCO with S < 0.15. 
|p^-p7| This is also the case in the recent measurements 
on the electron doped material NCCO. Although 
above Tc, the non- metallic behavior is clearly seen, be- 
low Tc, metallic-like behavior reminiscent of a Drude edge 
for wavevectors in the range 10 to 200 cm~^ has been ob- 
served, although none has yet been seen in Bi2212. 
p3[ This metallic-like behavior is most likely associated 
with the c-axis supercurrent. Thus, if Bi2212 were to 
behave the same as the other materials with an incoher- 
ent c-axis normal state conduction, one would expect it 
to show this same metallic-like behavior below Tc. 

There have been two schools of thought on this is- 
sue. One is that the interlayer tunneling processes change 
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dramatically from incoherent to coherent below Tc- p9| ] 
Our analysis of the c-axis twist experiments of Li et al. 
strongly contradicts this idea, since the quasiparticle tun- 
neling below Tc must be strongly incoherent. The second 
school holds that there is no quasiparticle tunneling be- 
low Tc due to an "orthogonality catastrophe", but that 
the only tunneling process occurs by the simultaneous 
tunnehng of pairs. |Q However, the fact that HgB2- 
intercalated Bi2212 has nearly the same Tc as does unin- 
tercalated Bi2212 argues strongly that this interlayer pair 
tunneling model for the occurrence of superconductivity 
is not correct. Q 

Thus, if the superconductivity arises from intralayer 
pairing, then the tunneling process must be tunneling by 
individual quasiparticles. The c-axis twist experiments 
provide strong evidence that this tunneling is incoherent, 
meaning that the momenta of the quasiparticles parallel 
to the layers before and after the tunneling process are 
not correlated, or random. This is in addition to any 
possible changes in the kinetic energy that might also 
occur. (H] 

With incoherent normal state c-axis tunneling, it is ex- 
tremely difficult to obtain a significant amount of c-axis 
critical current with an OP that averages to zero on the 
Fermi "surface" . Since there is a very strong consensus 
that the superconducting critical current in all HTSC 
is three-dimensional, incoherent tunneling in the c-axis 
direction necessarily implies that the OP must have an 
s-wave component, so that its Fermi "surface" average 
is non- vanishing at all T < Tc- This is also true in the 
organic layered superconductor K-(ET)2Cu[N(CN)2]Br, 
which has recently been shown to have extremely inco- 
herent c-axis propagation, both above and below Tc- | |3l| ] 
Combined with this supporting evidence that the c-axis 
tunneling in most HTSC is predominantly incoherent, 
the twist experiments of Li et al. provide very strong 
evidence of an s-wave superconducting OP in the HTSC. 

Very recently, a simple model for the apparent vio- 
lation of the conventional sum rule observed by Basov 
et al. was investigated theoretically by Kim and Car- 
botte (KC). 1^ Those workers assumed a dj.2_y2- 
wave OP, and an incoherent interlayer tunneling ma- 
trix element squared of the form /'^(0k,0k') = \Vq\ + 
iV^iP cos(2(/>k) cos(2(/>k')- By equating derived from 
the conductivity and the superfluid density ps , each eval- 
uated to lowest order in /•■', they wrote, [B2[ 



6N,{T/Tc)='- + ^-^^''^'^''''^'''^^^ 



]) 



2E.(/^('/'k,0k')'5ck^Ck' 



(13) 



where SNt^ = {N^ - iVsc)/PsC; ^Ck = ^^/^Ck, ^ck' = 

Cj/r2^k', ^Ck = Ai;((/)k)/f^Ck, (^Ck' = ((/)k')/f^Ck', ^^Ck = 

[^^2 + A2(0k)]^/^ andf^ck' = [w^.^ A^((/>k')]^''^- It is easy 
to show that for coherent tunneling SNq{T/Tc) — 1 for 
any OP, for all T/Tc < 1. It is also easy to show that 



the rf-wave case, KC showed that incoherent tunneling 
gives a conductivity sum rule violation that is strong and 
of the wrong sign. j2^j3^,Q They gave a lower limit on 
the violation, based upon restrictions of the parameters 
Vq and Vi. However, KC did not calculate the extended- 
s-wave case, which is actually very interesting to do. 

We studied the one-parameter model of Graf et al., 
pO| , which interpolates smoothly between coherent and 
incoherent tunneling. In this model, we can analytically 
perform the calculations at Tc for the = d,e cases. 
We find that dNd{l) = [2 + l//d(7)]/4 and SNdl) = 
[2 -I- l//e(7, 0)]/4. In the coherent limit 7 00, all three 
OP's give no sum rule violation, as expected. However, 
in the incoherent limit 7^0, SNj_{l) 4/7^, which 
diverges strongly, and SNc{l) ^ (1 + 7rV8)/2 w 1.117. 
As r ^ in the coherent limit 7 — > c», we again have 
SN(^{0) — 1 for ( ~ s,d, e, as expected. As T ^ in the 
incoherent limit 7 — + 0, for both C = d, e, we can evaluate 
the denominator in Eq. ( [l3| ) exactly, and the numerator 
numerically. We find ^iVe(O) 1.087 and SNd{0) 00. 
We thus conclude that none of the three OP's gives a 
SN < 1, but the d-wave case is by far the worst, espe- 
cially in the limit of strongly incoherent tunneling. Thus, 
it is likely to be much easier to construct a theory that 
can incorporate both incoherent interlayer tunneling and 
a SN < 1 if the OP is an anisotropic s-wave one, and not 
d-wave. In particular, theories based upon two gaps, a 
non-superconducting pseudogap and a superconducting 
gap, appear likely to satisfy the experimental sum rule 
violation, provided that the OP describing the supercon- 
ducting state is either s-wave or extended-s-wave. 

In conclusion, we found that the c-axis twist experi- 
ments of Li et al. provide compelling evidence that the 
c-axis tunneling in Bi2212 is strongly incoherent. Q| As a 
consequence, the experiment cannot distinguish between 
an isotropic s-wave and an extended-s-wave OP. How- 
ever, the purported c?2,2_j^2-wave OP can be ruled out in 
Bi2212. The other HTSC and organic layered supercon- 
ductors which also have such incoherent c-axis tunneling 
also cannot have a dx2-y2-'wave OP. Such incoherent c- 
axis tunneling is likely to be consistent with the available 
c-axis infrared reflectivity measurements, provided that 
the OP is not of the dj.2_j^2-wave form. 
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nanostrukturierter Festkorper." 



5N^ (T/Tc) = 1 , regardless of the form of Z'' ((^k , (j^k' 
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